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By analyzing about 33 pb^^ data sample collected at and around 3.773 GeV with the BES- 
II detector at the BEPC collider, we directly measure the branching fractions for the neutral 
and charged D inclusive semimuonic decays to be BF{D^ l^^ X) = (6.8 ± 1.5 ± 0.7)% and 
BF{D^ — > X) = (17.6 ± 2.7 ± 1.8)%, and determine the ratio of the two branching fractions to 
be |^gr:+^)' = 2.59 ±0.70 ±0.25. 
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I. INTRODUCTION 

The neutral and charged D arc both charmed mesons. 
However, the hfetime of the meson is surprised longer 
than that of the meson [l| . Isospin symmetry predicts 
that the partial widths to Cabibbo-favored semileptonic 
decays of the and mesons are equal. It is expected 
that the ratio of the branching fractions for the D'^ and 
D"*" inclusive semileptonic decays is approximately equal 
to the ratio of D° and D'^ lifetimes up to the order of 
O{0^) Measurements of the branching fractions for 
the D° and D'^ semileptonic decays can provide valu- 
able information on the difference in their lifetimes, which 
is important to understand the charmed meson decays. 
Moreover, by comparing the branching fractions for the 
and inclusive decays with the sum of those for 
the exclusive decays, one can estimate whether there are 
some decay modes which are not observed yet. The D 
mesons inclusive semielectronic decays have been studied 
by several experiments However, the experimen- 

tal studies for the D mesons inclusive semimuonic decays 
are very limited. Actually, there is no measurement for 
the decay £)+ — > fi^X available (X=any particles) in the 
PDG 111 now. 

In this Letter, we report direct measurements of the ab- 
solute branching fractions for the inclusive decays — > 

X and fx'^X , based on analyzing a data sample 

of about 33 pb~^ collected at and around the center-of- 
mass energy (^s) 3.773 GcV with the BES-II detector. 



the total solid angle with longitudinal (transverse) spa- 
tial resolution of 5 cm (3 cm) . End-cap time-of-flight and 
shower counters extend coverage to the forward and back- 
ward regions. A Monte Carlo package based on GEANT3 
has been developed for BESII detector simulation and 
comparisons with data show that the simulation is gen- 
erally satisfactory 0] . 

III. DATA ANALYSIS 

Around 3.773 GeV, the ■i/'(3770) resonance is pro- 
duced in electron-positron (e~^e~) annihilation and de- 
cays into DD pairs with a large branching fraction, 
BF(V'(3770) ^ DD) = (85 ± 5)% Q, i, i, [Ol. If 
we can reconstruct D (they are named as singly tagged 
D) from the DD pairs, the other D must exist on the re- 
coil side of the tagged D. Taking this advantage, we can 
directly measure the absolute branching fractions for the 
D decays with the singly tagged D samples. The singly 
tagged D samples used in the analysis were reconstructed 
from the hadronic decay modes of mKnir (m=0,l, 2; 
ri=0, 1, 2, 3, 4) in the previous works [l^, These 
give the total numbers of 7584 ± 198(stat.) ± 341(sys.) 
singly tagged D^ [ll and 5321 ± 149(stat.) ± 160(sys.) 
singly tagged D~ ^M,- Throughout the Letter, charge 
conjugation is implied. 



A. Candidates for D ^ X 



II. THE BEIJING SPECTROMETER 

BESII is the upgraded version of the BES detec- 
tor operated at the Beijing Electron Positron Col- 
lider (BEPC) 0. A 12-layer vertex chamber (VC) sur- 
rounding the beam pipe provides trigger information. 
A forty-layer main drift chamber (MDC) located out- 
side the VC performs trajectory and ionization energy 
loss (dE / dx) measurement with a solid angle coverage 
of 85% of in for charged tracks. Momentum resolution 
of (Jp/p = 1.7%-\/l + p^ {p in GeV/c) and dE/dx resolu- 
tion of 8.5% for Bhabha scattering electrons are obtained 
for the data taken at ^/s = 3.773 GeV. An array of 48 
scintillation counters surrounds the MDC and measures 
the time of flight (TOF) of charged tracks with a res- 
olution of about 200 ps for the electrons. Surrounding 
the TOF is a 12-radiation-length, lead-gas barrel shower 
counter (BSC) operated in limited streamer mode, which 
measures the energies of electrons and photons over 80% 
of the total solid angle, and has an energy resolution 
of cfe/E = 0.22/v^ {E in GeV), spatial resolutions of 
(T0 = 7.9 mrad and az = 2.3 cm for the electrons. Out- 
side of the BSC is a solenoidal magnet which provides a 
0.4 T magnetic field in the central tracking region of the 
detector. Three double-layer muon counters instrument 
the magnet flux return, and serve to identify muon with 
momentum greater than 0.5 GeV/c. They cover 68% of 



The candidates for D /i+X are selected in the sys- 
tem recoiling against the singly tagged D mesons. It is 
required that the candidate tracks should be well recon- 
structed in the MDC with good helix fits, and satisfy 
|cos0| < 0.67, where Q is the polar angle. Each track 
must originate from the interaction region, which requires 
that the closest approach to the interaction point in the 
xy-plane is less than 2.0 cm and in the z direction is 
less than 20.0 cm. To reject the muons from kaon and 
pion decays, the selected candidate tracks are required 
to originate from the same vertex as those decay from 
the singly tagged D mesons, by requiring Sz < 2az (2.0 
cm), where dz is the minimum distance in the z direc- 
tion between the candidate track and those decay from 
the singly tagged D, az is the standard deviation of the 
Sz distribution. To ensure that the track can be detected 
by the muon counter, the transverse momentum of each 
track is required to be greater than 0.52 GcV/c. Since 
the hit depth of muon in the muon counter varies with its 
transverse momentum, the muon is selected by requiring 
that the track should hit at least one layer for the track 
with the transverse momentum in the region (0.52, 0.75) 
GeV/c, two layers for the one in the region (0.75, 0.95) 
GeV/c, or three layers for the track with transverse mo- 
mentum being greater than 0.95 GeV/c, respectively. Be- 
cause there is no charge symmetric background as those 
in the measurements of the semielectronic decay 0, 0], 
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only " right-sign" muon candidates (with their eharge op- 
posite to the flavor of the single tag) are selected. 

Fig. [1] shows the resulting invariant mass spectra of 
the mKmr combinations for the events in which the can- 
didates for muon are observed on the recoil side of the 
mKn-K combinations for the singly tagged (left col- 
umn) and (right column) mesons. Fitting each in- 
variant mass spectrum with a Gaussian function for D 
signal and a special function [13] to describe the back- 
ground shape, we obtain number A^^j^^ of the observed 
muon candidates. The total numbers of the observed can- 
didates for D° ^i+X and £>+ ^ ^l+ X are 79.3 ± 10.3 
and 99.6 ± 11.7, respectively. 
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FIG. 1: Invariant mass spectra of the mKnn combinations 
for the events in which the muon candidates are observed on 
the recoil side against the mKmr combinations for the singly 
tagged (left column) and (right column) mesons, 
where the tracks are with the transverse momenta in the re- 
gion: (a) and (f) (0.52, 0.62) GeV/c; (b) and (g) (0.62, 0.72) 
GeV/c; (c) and (h) (0.72, 0.82) GeV/c; (d) and (i) (0.82, 0.92) 
GeV/c and (e) and (j) (0.92, 1.02) GeV/c intervals, respec- 
tively. 



B. Unfolding Procedure 

Since the muon, electron, kaon and pion may be 
misidentified to each other, the observed muon candi- 
dates consist of the true muons and those are misidenti- 
fied from electron, kaon and pion. The true yield N^,.^^ of 
muons can be extracted through an unfolding procedure. 
A detailed description of the unfolding procedure can be 
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found in Ref. [J • The right-sign candidates are unfolded 
using the matrix Eq. [l] 



true 
^ *tru( 

(1) 

where N^-^^ and N^^^^ [a and h denote muon, electron, 
kaon or pion) are the number of the observed candidates 
for a and the number of the true a particles, respectively; 
^b^a YsXe of misidentifying the particle as a (a ^ 

h) or the efficiency of identifying the particle a {a = b). 

To remove the misidentified particles from the observed 
muon candidate sample, we need to know the numbers 
A^Qi-,g, and N^^^^ of electrons, kaons and pions be- 

sides the number iV^j^g of muons. They are also se- 
lected in the system recoiling against the singly tagged 
D mesons. Electron, kaon and pion are identified by us- 
ing the dE/dx, TOF and BSC measurements, with which 
the combined confidence levels for the electron, kaon and 
pion hypotheses {CL^, CLk and CL^r) are calculated. 
An electron candidate is required to have CL^ > 0.1% 
and CLe/ {CLe + CLk + CL.^) > 0.8, and a kaon candi- 
date is required to satisfy CLk > CL.^ and CLk > 0.1%. 
The tracks not satisfying the selection criteria of muon, 
electron and kaon are treated as pions. 

To obtain the number of the true muons, we also need 
to estimate the rates ry''"*" by analyzing pure muon, 
electron, kaon and pion samples. The muon sample is 
selected from cosmic rays. The electron sample is se- 
lected from the radiative Bhabha events. The pion and 
kaon samples are selected from the J/ip unr'^TT~ and 
J/t/j (j>K~^K~ events, respectively. The difference be- 
tween the DD decay event environment and the selected 
particle sample environment is studied with two Monte 
Carlo matrices. One matrix is determined using the 
muons. electrons, kaons and pions from the Monte Carlo 
simulated cosmic ray, Bhabha. J/ip (f>K'^K~ and 
J/ip ^ LL)Tr^TT~ decays, and another matrix is produced 
using the particles from the DD Monte Carlo samples. 
The difference in the unfold yields with various matrices 
are about 3.9% for D^ decay and 2.1% for decay. 
Since the rates ry*^" vary with momenta of the particles, 
we divide the transverse momentum region (0.52,1.02) 
GeV/c into five intervals in the analysis. Figure shows 
the rates each transverse momentum interval ob- 

tained from the pure particle samples. 

Inserting the numbers of iV",^g, the rates rf'~*°->^ {i de- 
notes the zth transverse momentum interval) in the ma- 
trix Eq. [Tl we obtain the yield N^^ 
the zth momentum interval. Summing the N^^\ 
intervals, gives the total numbers of the true yield candi- 
dates to be 87.5 ± 17.1 for £)° X and 129.3 ± 19.5 
for _D+ M^AT, respectively. 



/^.^Q of the true muons in 
over the 
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IV. RESULTS 
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FIG. 2: rj''^" {a and b denote muon, electron, kaon or pion) 
are the rates of misidentifying the particle b as a {a b) 
or the efficiency of identifying the particle a {a = b), where 
the dots, squares, triangles and stars are for the particles 
(mis)identified as muon, electron, kaon and pion from the 
pure particle samples: (a) the muon sample selected from the 
cosmic rays; (b) the electron sample; (c) the kaon sample and 
(d) the pion sample. 



by 



The branching fraction for D ^~^X is determined 



BF{D n+X) = 



Nr 



ii+X 



(2) 



where Njj^^+x is the number of the events for D 
fi^X, N^, is the number of the singly tagged D mesons, 
Cd^^i+x is the detection efficiency for D /z+X. 

The detection efficiency e^^^+x is estimated with the 
Monte Carlo simulation. The Monte Carlo events are 
generated as e+e~ DD, where D decay into the singly 
tagged D modes and D decay into semimuonic modes. 
The semileptonic decays arc generated with the de- 
pendence of form factors given by the pole model [l^ . 
The generator has been applied in the previous mea- 
surements and reproduces exclusive semileptonic spectra 
well [H, [H, [13, til, [il- In order to study the model 
dependence of the efficiency, the semileptonic decays are 
also simulated with the ISGW2 form factor model [20| . 
The difference in efficiencies determined with various 
form factor models is about 5.8%. The averaged effi- 
ciency is determined by weighting the branching frac- 
tions of the D meson semimuonic decays [i| and the 
numbers of the singly tagged D events. The efficiencies 
are (15.4 ± 0.2)% for D° and (13.5 ± 0.2)% for 

fi~^X, respectively. 

Inserting the number Njj^^+x of the signal events 
for D l^'^X, the number Njj of the singly tagged D 
mesons and the detection efficiency e^^^+x in Eq- ©j 
we obtain the branching fractions for D ^ ji^X iohe 



BF{D" ^i+X) = (6.8 ± 1.5 ± 0.7)% 



and 



C. Other Background 



In the selected candidate events, there are still some 
contaminations from decays M^^^m ^'^'^ ^ 

fj.'^v^. These backgrounds are estimated based on the 
Monte Carlo simulation. The simulated background 
events are generated as e+e~ DD events, where D 
and D mesons are set to decay into all possible modes 
except D'^ fi^X for studying the 1?° decay, or ex- 
cept D^ — > ji^X for studying the Z?+ decay, respec- 
tively. In the Monte Carlo simulation, the decay modes 
of D mesons and their branching fractions are quoted 
from PDG [l|, and the particle trajectories are simulated 
with the GEANT3 based Monte Carlo simulation pack- 
age of the BES-II detector 0. The numbers of these 
background events are estimated to be 8.4 ± 1.8 and 
3.0 ±1.6 for D" fi+X and D+ ^ n+X, respectively. 
After subtracting these background events, we obtain 
79.1 ± 17.2 and 126.3 ± 19.6 signal events for D° n+X 
and fJ'^X, respectively. 



BF{D+ ii+X) = (17.6 ± 2.7 ± 1, 



where the first errors are statistical and the second sys- 
tematic. 

The systematic errors arise mainly from the uncer- 
tainties in particle identification (~ 5% for muon [2l|). 
in tracking (~ 2.0% per track), in the number of the 
singly tagged D mesons (~ 4.5% for 5° [H and - 3.0% 
for D~ [IJl), in the Sz selection criterion (~ 3.5%), in 
the input form factor models 5.8%), in the selected 
sample environment (~ 3.9% for and 2.1% for D^), 
and in the Monte Carlo sample statistics (~ 1.5%). The 
systematic errors from the uncertainties of the 
the unfolding procedure are estimated using the Monte 
Carlo samples generated with the Gaussian distribution 
to describe the uncertainties of the r]''^°- matrices, and 
they are estimated to be ~ 2.1% for D° and ~ 1.0% for 
-D"^, respectively. The systematic errors from the poorly 
measured exclusive semileptonic decay modes are esti- 
mated using DD Monte Carlo samples generated with 
or without those modes, they are about 3.3% and 4.7% 
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for D° l^^v^ and £>+ M^^M' respectively. Adding 
these uncertainties in quadrature yields the total system- 
atic errors to be 11.0% and 10.3% for £)° ^i+X and 
D'^ fj,~^X, respectively. 

With the measured branching fractions for Z?" ji^X 
and D'^ — > n'^X, the ratio of the two branching fractions 
is determined to be 



BF{D+ -> fi+X) 
BF{D° li+X) 



2.59 ±0.70 ±0.25, 



where the first error is statistical and the second system- 
atic arising from the uncertainties in the number of the 
singly tagged D mesons, in the rf'^"-^ in the Monte Carlo 
sample statistics, in the the selected sample environment 
and in the poorly measured decay modes. 

TableUshows the comparisons of the measured branch- 
ing fractions for D yU+X by the BES Collaboration 
with those measured by the ARGUS [13, CHORUS il 
Collaborations and the averaged value from PDG [l|. 
The measured BF{D^ IJ-^X) is in good agreement 
with the measurements from other Collaborations. 



V. SUMMARY 

Using the data sample of about 33 pb~^ collected at 
and around ^/s = 3.773 GeV with the BES-II detec- 



tor at the BEPC collider, we have studied the inclusive 
semimuonic decays of D mesons. The absolute branch- 
ing fractions for /i+AT and fi~^X are mea- 
sured to be BF{D" fi+X) = (6.8 ± 1.5 ± 0.7)% and 
BF{D+ fi+X) = (17.6 ± 2.7 ± 1.8)%. The latter one 
is the first measurement. With the measured branch- 
ing fractions for 1?° /i+A" and fJ-~^X, the ra- 
tio of the two branching fractions is determined to be 
BF{D+ ^i+X)/BF{D° ^i+X) = 2.59 ± 0.70 ± 0.25, 
which is consistent with the ratio of the lifetimes of D'^ 
and D° mesons, r£,+ /T£,n = 2.54 ±0.02 Q. 
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TABLE I: Comparisons of the measured branching fractions for the inclusive semimuonic decays of D mesons with those 
measured by ARGUS [H, CHORUS [H Collaborations and the PDG values f^. 



BES ARGUS CHORUS PDG06 

BFiD" tJ.+X)[%] 6.8±1.5±0.7 6.0±0.7±1.2 6.5±1.2±0.3 6.5±0.7 
BF{D+ ^ lj,+X)[%] 17.6±2.7±1.8 

B^ffr:+g 2.59±0.70±0.25 

^ - - - 2.54±0.02 



